Preprint typeset in JHEP style. - PAPER VERSION 



Bicocca-FT-99-33 
November 1999 



Power corrections to the differential 
Drell-Yan cross section 



Mrinal Dasgupta 

Dipartimento di Fisica, Universitd di Milano Bicocca 

and INFN, Sezione di Milano, 

Via Celoria 16, 1-20133, Milano, Italy 

E-mail: dasguptaSmi . inf n . it 



Abstract: We estimate the power corrections (infrared renormalon contributions) 
to the coefficient functions for the differential Drell-Yan cross-section (Pa/dQ'^dy, 
where is the mass squared and y the rapidity of the produced lepton pair. We 
employ the dispersive method based on the analysis of one-loop Feynman graphs 
containing a massive gluon. 
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1. Introduction 

The Drell-Yan process [1] describes the colhsion of two hadrons and the subsequent 
production of a lepton pair and a hadronic final state. In the parton model it proceeds 
simply via the annihilation mechanism where a quark and anti-quark generated by 
the parent hadrons annihilate to form a photon which decays to a lepton pair. At 
present day collider energies one is also able to produce the electroweak W and Z 
bosons on mass shell via this mechanism. 

Historically the Drcll-Yan process has played an important role in the develop- 
ment of QCD. The parton model became more firmly established when it was realised 
that it gave a good description of the data in hadron-hadron collisions. In addition 
it became evident that QCD perturbation theory could be applied to describe strong 
interaction phenomena when one encountered the very same mass singularities in 
Drell-Yan calculations as those in the case of deeply inelastic scattering giving rise 
to the concept of universal functions that control the long-distance dynamics, which 
one refers to as parton distributions. 

Data from hadron-hadron collisions has proved a valuable source from which one 
has been able to constrain and measure various parton distributions. In particular 
such data provides a means to extract information on the quark distributions in 
pions which is inaccessible from DIS experiments. Data on low-mass lepton pair 
production has been used to study the small x behaviour of parton distributions. 
For a thorough review of available Drell-Yan data and comparisons with theory the 
reader is referred to Ref . [2] . 

In this article we choose to concentrate our attention on the cross section dPa/dQ^dy 
with Q'^ being the mass squared and y the rapidity of the produced lepton pair. Most 
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recently the CDF collaboration have studied the rapidity dependence of the Drell- 
Yan cross section in a limited rapidity range [3]. In fact from the experimental side 
there is a wealth of data on various rapidity distributions (see Ref. [2] ) but progress 
on the theoretical side is somewhat lacking. While the 0{as) perturbative QCD 
calculations for the above distribution were performed several years ago [4, 5] there 
is as yet no 0(0^) estimate. Also lacking is any estimate of the non-perturbative 
power-like corrections which have been extensively studied in many other cases . The 
aim of this article is to study the power correction to the rapidity distribution which 
should in principle allow a better description of the data when added to the pertur- 
bative predictions. Power correction predictions already exist for the more inclusive 
cross section da/dQ^ where a l/Q^ dependence is predicted with a characteristic 
phase space enhancement [6,7] . 

Prom a purely theoretical viewpoint the motivation of the work described here is 
the testing of current ideas on power corrections which seem to indicate that though 
these contributions are non-perturbative in origin one may suitably extend a pertur- 
bative approach to estimate them. The success of such an approach in the case of 
DIS structure functions [8] is encouraging enough to extend this study to the present 
case. In view of the relative simplicity of the renormalon calculations that yield 
predictions for the power corrections and the range and accuracy of current exper- 
imental data on several different QCD observables one should be able to undertake 
a serious and extensive confrontation of these theoretical ideas with the data, a task 
that has already begun [9-11]. 

This paper is organised as follows. In the next section we give a very brief 
review of the dispersive treatment of power corrections which has been described in 
great detail previously [6] . Following this we mention the notation and introduce the 
kinematical variables relevant to our study. In the following section we describe our 
results for the power corrections and finally make some concluding remarks. 

2. Dispersive Approach 

The main ideas of the dispersive approach to power corrections can be briefly summed 
up as below. First one assumes a QED inspired dispersion relation to be formally 
true in the QCD case so that 

"^(^') = -r^rpp^(/^') (2-1) 

Jo jl'^ -\- /c^ 

with the spectral function 

Ps{l^') = ^ {asil^'en - asil^' e— )} . (2.2) 

Thus one is assuming that the QCD coupling is well behaved in the infra-red and 
the only singularity is a discontinuity on the negative real axis of its argument. 
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Non-perturbative effects at long distances are expected to give rise to a mod- 
ification to the perturbatively-calculated strong coupling at low scales, 6as{fJ.'^) = 



a; 



a; 



(yU^) being the perturbatively-calculated running coupling [6]. 



The spectral function (2.2) receives the corresponding modification 



(2.3) 



To consider the effect of the above on an observable F one assumes the implicit 
inclusion of a gauge invariant set of higher order graphs (which in a large Nf ap- 
proximation just reduce to quark bubble insertions ) combined with single gluon 
emission has the effect of generating the running coupling in the one-loop calculation 
of -F . In practice the running coupling can only be reconstructed in this manner, 
for a sufficiently inclusive observable like the one we study in this article. The above 
considerations allow one to write [6] 



(2.4) 



where JF is the one loop correction to the observable, computed with a finite gluon 
mass \j? . Then the non-perturbative region contributes to F through Eq. (2.3) and 
one finds 



?>F 



'7 



fc,(/x2)^(/x7g2) 



where setting /x^/ Q"^ — e 
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6;(e) = -— Disc{^(-e)} . 



(2.5) 



(2.6) 



Since 6as{fJ^'^) is limited to low values of /i^ the asymptotic behaviour of 6F at high 
is given by its behaviour in the limit e — 0. Clearly only terms that are non- 
analytic in e in the small e behaviour of J^, yield non-perturbative modifications to 
F within the above approach. In the present case the characteristic function will 
be found to have the general small e behaviour 



^ ^IT- (CiMeln^e + CsMelne) 
27r ^ ' 



(2.7) 



with {x} denoting the phase space dependence that we compute here. According to 
(2.5) and (2.6) the above e dependence translates into the power correction 



C, {x} + 2Ci {x} In f ^ 



with the parameters ^2 and B2 being 

Cf d^x^ 
2^ 



A2 



Jo 

J_Cf r 

A2 27r Jo 



1^' 



(2.8) 



(2.9) 



fj,^ In 11^ 5q;s(a*^) • 
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Then one invokes the universahty assumption, in that if one beUeves that the concept 
of as can be meaningfully extended to small scales, it can be done in an observable 
and indeed process independent fashion. This would allow us to extract the value of 
the above defined moments of the coupling, A2 and B2, from any experimental data 
where the same power correction is obtained and use it to fit Drell-Yan data. Recent 
studies of 1/Q corrections to event shape variables in e~^e~ annihilation [10] and in 
DIS [11] lend support to this idea . Studies of the non-singlet contribution to 

DIS structure functions suggest that A2 — 0.2 GcV^ [8]. We do not know the value 
of the parameter B2 from any experimental data and hence it can be treated as a 
free parameter in the fit to the Drell-Yan case. 

On the other hand as has been explained in detail in previous articles (see for 
instance Ref. [12]) the parameter A2 and indeed any other moment of the modified 
coupling, 5 as, depends on the order of the perturbative result with which one is 
interested in merging the non-perturbative piece. The above DIS value is relevant to 
merging with the next-to-leading order (NLO) QCD prediction while in the present 
case of the differential Drell-Yan cross-section one only knows the leading QCD cor- 
rection to the parton model. However we use the DIS value for illustrative purposes 
here and in the expectation that the NLO result will be available soon. 



3. Definitions and Kinematics 

We wish to compute the power corrections to the observable (P a / dQ"^ dy where is 
the mass squared and y is the rapidity of the produced lepton pair. To lowest order 
the expression for the above quantity is simply 



dQHy 



Ns 



(3.1) 



where 



N is the number of colours and s the center-of-mass energy squared. In addition 
Cg refers to the charge of fiavour q in the fiavour sum above and Zi, Z2 refer to the 
momentum fractions carried by the quark and antiquark, of the parent hadrons A 
and B while the / functions are the corresponding parton distributions. We explicitly 
have, in the hadronic centre of mass frame, the parton momenta (neglecting their 
small intrinsic kt) 



fs 

Pi = ^(^1,0,0,^1) 

P2 = —(2:2,0,0,-^2) 



(3.2) 
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The invariant mass of the produced lepton pair is then simply given by 



<5 =(Pl+P2) =SZi^2 (3.3) 



'.^(^)-h^(^) (3.4) 



while its rapidity is 
which yield 

Zi — Vre^ , Z2 — Vre ^ , r = — . (3.5) 

Beyond the naive parton model one has to consider radiative corrections to the above 
picture, wherein initial partons carrying momentum fractions zi and Z2 degrade their 
momenta via gluon radiation before annihilating to form the lepton pair. Hence 
the simple correspondence (3.5) between the momentum fractions of the initiating 
partons and the mass and rapidity of the lepton pair no longer holds beyond lowest 
order. We define for later use kinematic variables Xi and X2 such that 

xi = Vre^ , X2 = \fTe"^ (3.6) 

and which at Born level arc just equal to the initial parton momentum fractions. 

For computing the rcnormalon contribution to the desired observable we shall 
need the squared matrix element for the QCD radiative process ^ 

?(pi)+g(P2)^7*(g)+^(^), (3.7) 

which takes the simple form 

where s,t,u are Mandelstam invariants defined as 

s = (Pi +P2)', i^{pi - k)\ u^ipi- qf, (3.9) 
and which satisfy the relation 

s + i + u^Q^ + l_i'^ (3.10) 

with /i^ = eQ^ being the squared gluon mass and s = Z1Z2S. 

The rapidity y can readily be obtained in terms of the Mandelstam invariants 
and the momentum fractions Zi and Z2 as 

y^UJ (3.11) 
^ 2 \s + t- Z2) ^ ' 



^For computing the power correction we are only considering tlie anniliilation contribution. In 
general there is also the QCD Compton scattering process to consider, beyond the naive parton 
model. We shall comment further on this in the final section. 
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which is equivalent to 

i = Q^^2^i -xiZ2{s- n'^) 

X2Z1 + X1Z2 

u — 

X2ZI + XIZ2 

with xi and X2 defined as before. The phase space for massive gluon emission is 
given by 

{zi - Xi){z2 - X2) > — . (3.13) 
Putting the above together we obtain the required differential cross section 

(Pa T{T + ZiZ2-eT) , , 

,^2^ TT , M{s,Xi,X2,Zi,Z2,€) (3.14) 

dQ^dy {ziZ2)[ziX2 + 2:2X1)^ 

where M. is just the matrix element (3.8) with the substitutions of Eq. (3.12) and A — 
16a'^ a sCg/ 27 Q'^s. The hadron level result is related to the corresponding partonic 
quantity by folding with parton distribution functions as below 



dQ^dy ?ll/"^^"^dg4®r^""^^^'^""^^"^j^^^'^''^^ ^'-''^ 
where for brevity we defined 

:Fg{z^,Z2) = {fg/A{zi,Q')f^/B{z2,Q') + {q^ q)} • (3.16) 

Next we introduce the variables ^ = xi/zi and ( — X2/Z2 in terms of which 
Eq. (3.15) assumes the familiar form 

-Y'^^ C ficiUMiM. - 0(1 - C) - <C) 

(3.17) 

in writing which wc have dropped an uninteresting overall constant factor a^/Ns, 
which also appears in the Born cross section to which we shall normalise the result 
subsequently. The coefficient function C(^, C, e) then takes the form 

r(^^A 2(1 +ec- 6^0(1 +e^c^(i + e)^) , ec(i+ec-6ec) 



2eeC(l + eC-eeC) 



e e 
+ 



Next taking moments of the convolution equation (3.17) we find as usual (note 
the different powers) 

dx, dx2 x^'^^ = C{N, M, e) ^,{N, M) (3.19) 
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where 

C(N, M, e) = f f C(e, C, e) e((l - ^(l - C) - e^Qdi dC (3.20) 

JO JO 

T,{N, M) = f,/A{N) fq/sm + (g ^ q) 

and the / functions represent the MeUin transforms of the quark and anti-quark 
density functions. 

As explained previously, to extract the power corrections we have to look for non- 
analytic behaviour in e, in the small e expansion of the Mellin transformed coefficient 
functions C{N,M,e). This non-analyticity will manifest itself (in the present case ) 
through the appearance of a logarithmic dependence on e in addition to the usual 
logarithmic divergences generated by the soft and coUinear regions of integration. 
These divergences are of course cancelled by virtual corrections in the case of the 
infrared and by absorbing the coUinear divergences into the definition of the parton 
densities. We are then left with terms like e In^ e and e In e which are in one to one 
correspondence with infrared renormalon poles in the Borel plane and will generate 
power corrections. We neglect similar terms that appear at order and higher 
orders in e since they will induce sub-leading 0{1/Q'^) power corrections in which 
we are not phenomenologically interested. 



4. Power Corrections 

In taking the Mellin transforms of the coefficient functions as mentioned in (3.19) 
one finds that the second and third pieces of the expression on the right-hand-side 
of Eq. (3.18) do not produce any logarithmic divergences but contribute only to 
the power corrections through the appearance of an elne term. The coUinear and 
infrared divergences lie in the first piece on the right-hand-side of (3.18). The Mellin 
transforms of this piece are cumbersome to evaluate directly and the result is most 
easily arrived at after a further change of integration variables: 

r = ^C,V = ^- (4.1) 

In terms of these variables the double MeUin transform of the above mentioned 
piece takes the simple form 

r''"' dt r f '"^""^ (4.2) 

JO Jr;i {1 — riT — €T){ri — T — €r]T) 

In the above equation the limits of the rj integration are given by 



r/i = — = — (A + 2f + VA2 + 4Af ), A = 1 - 2f + - 2fe - 2f^e + e^f^. (4.3) 
ri2 2t 
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Taking the 1] moments is now an easy task but inserting the hmits (4.3) comph- 
cates the extraction of the r moments. In particular one has to evaluate the form 
(apart from other relatively straightforward integrals) 



i/(i+VJ)' 



dr 



1 + t' + f'e 



2^2 



2ef2 



1 -f 



tanh 



-1 



er 



1 — f — ef 



(4.4) 



where we use uj to denote a generic power which depends on M and N. The above 
form has been evaluated in Ref. [6] and we shall not describe the manipulations that 
lead to the evaluation of integrals of the above type but refer the reader to section 
4.6 and appendix A of Ref. [6] for the details. The moments of the second and third 
pieces on the right hand side of Eq. (3.18) are most easily evaluated directly in ^, ^ 
space. 

Putting together the contribution of all the pieces and including the time-like vir- 
tual corrections (which arc independent of M and N and were computed in Ref. [6]) 
we obtain the result for the total contribution as follows: 



C^+'^iN, M, e) = Co In e + e In^ e + In e 



where 



Co 

Ci^- 



2Si{N + 1) + 25i(M + 1) - 3 -h 



1 



1 



1 



1 



Co 



M+1 N+1 N + 2 M + 2 

N M 

[-{M + N + 2)Si{N +1)-{M + N + 2)Si{M + 1) + 2MN + bM + bN 
N-M M-N N-2M M-2N 



+ 



+ 



+ 1 M+1 
- 2M - 3 M 

+ ^TTTT ^ + 



+ 



+ 



2(A^ + 2) 2(M + 2) 
2A^-3" 



(4.5) 

(4.6) 

(4.7) 
4 
(4.8) 



2(A^ + 3) 



2(M + 3) 



and 



N-l 



Si{N)= Y,^ = ^{N)+^E. 
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(4.9) 



Notice the absence of a In^e term which cancels in the sum of the real and 
virtual pieces. In addition, the logarithmic divergence in the cut-off (gluon mass) is 
absorbed into the structure functions and will not concern us any more. We shall 
instead concentrate our attention on the corrections generated by the e In^ e and e In e 
piece of Eq. (4.5). As a check on our result above, we note that in the diagonal limit 
M = N we recover the result obtained in Ref. [6] provided one changes A^ to A^ — 1 
in accordance with the different definition of the moments adopted in that reference. 
This agreement is expected as in the special case, M — N, our calculation reduces 
to just extracting the r dependence of the power corrections. The above results 
can be easily expressed in ^, ^ space and the convolutions with the parton densities 
performed as prescribed in Eq. (3.17). 
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5. Discussion 

For the purposes of illustration we assume that fg/A — fq/B — Q with the q q term 
being similarly labelled q. In reality of course such assumptions about the density 
functions will depend on the beam and the target used in the relevant experiment 
but our predictions can be easily adjusted to every case. 

600 



500 
400 
5 300 
200 
100 





-3-2-10123 

y 

Figure 1: Plot of the rapidity dependence of the power correction function if(r, y). The 
solid line corresponds to r = 0.0149, |y| < 2 while the dashed line corresponds to r = 
0.0019, \y\ < 3. 

Then we find that performing the convolution and extracting the lowest (Born) 
order result allows one to write : 

(_fo_\ 
dQHy~\dQHy)^ 

where the power correction function H{xi,X2) can be expressed as 
H{xi, X2) = -rr-^ r ^qi {^1 + X2 q'{x2)q{xi) - 4 q{xi)q{x2)} In 

J[Xi,X2) q 

+ 2xiX2 q'{xi)q'{x2) - Sxi q'{xi)q{x2) - q'{x2)q{xi) 

- q{xi) I{x2) - q{x2) I{xi) - Xi q'{xi) K{x2) 

- X2 q'{x2) K{xi) + q'^ q] 

(5.2) 
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1 + ^ H{Xi,X2) 



(5.1) 




with the following definitions: 



f{xi, X2) = eg(?(2^i)?(2^2) + 



(5.3) 




) 



(5.4) 



(5.5) 



In the above formulae the '+' distributions have their usual meaning, the prime 
symbol on the q denotes a derivative and the scale (Q^) dependence of the parton 
distributions is understood. Note that the symmetry under the exchange of Xi and 
X2 implies that the power correction is a symmetric function about y — with our 
assumption about the parton densities, which would be valid in the pp case. For a 
general beam and target, asymmetry will be induced purely by the differing parton 
densities in the beam and target. 

A plot of the function H against the rapidity y is shown in figure 1, for two 
different Q values which correspond to two different values of r. The value of ^/s 
was chosen to be 1.8 TeV, and the plots were made using the MRSA valence parton 
distributions [13]. The value of the parameter B2 was chosen to be 1 GeV^ for 
illustrative purposes; in principle there is no reason why it should have a value close 
to the one chosen here. 

Both curves shown reflect a similar behaviour, namely that the power correction 
is quite flat until one starts to reach the edge of the rapidity range shown in either 
case. For example, in the plot at 80 GeV as one gets closer to y = 3, we start 
approaching the region where Xi = ^fre^ is near unity. In fact for y = 2.9 and 
Q = 80 GeV one finds Xi = 0.80. Similarly as one progresses towards more negative 
rapidity values, X2 = \fTe^'^ starts approaching unity. When either Xi or X2 gets close 
to unity an explosive behaviour of the power correction is witnessed. This behaviour 
is a reflection of the singular nature (as ^, C — ^ 1) of the derivatives of delta functions, 
which are obtained on inverse transforming the Mellin-space coefficient functions C\ 



At the edge of the rapidity range for the 80 GeV plot, the overall effect of the 
power correction A2H(T^y)/Q'^ is large enough to be comparable to ^^(Q) which 
makes it an important effect to consider while comparing perturbative predictions 
with the data. For larger Q values, although there is a logarithmic enhancement of 
the power correction in this case, the suppression by should reduce the significance 
of the correction. 

In the more inclusive case of the cross-section differential in the explosion 
of the power correction will only be important as both X\ and X2 approach unity, in 
other words in the limit r — > 1 [6] . This region is probably beyond any experimental 
interest as the Q value is too high (close to the centre of mass energy) in that case. 



and C2. 
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Hence while the power correction should not be an important consideration in the 
distribution, its presence should be felt in the combined Q^, y distribution, especially 
at moderate Q values and towards the edge of the allowed rapidity range. 

Lastly we comment on the fact that at higher rapidities one would expect the 
QCD Compton scattering process to become significant. We have not taken this 
into account here as it is not yet completely clear how to treat renormalon contribu- 
tions from incoming gluons as a genuinely gluonic contribution. The only treatment 
suggested till now is to compute the renormalon contribution, treating the gluon as 
an internal line radiated off an incoming quark which has been done for singlet DIS 
contributions (sec Rcf. [14] and references therein). This procedure can in principle 
be applied here in the Drell-Yan case. However it leaves the question of how one may 
unambiguously factor off the quark to gluon splitting in order to be able to convolute 
with the gluon density, which rapidly grows at small momentum fractions x. Till 
this issue is better understood we postpone further discussion on this topic. 
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